Interference and interaction effects in adiabatic pumping through quantum dots 
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In order to investigate the effects of interference and interaction in adiabatic pumping, we consider 
an Aharonov-Bohm (AB) interferometer with a quantum dot embedded either in one or in both arms. 
We employ a real-time formalism and we perform an expansion both in the tunnel-coupling strengths 
between dot and leads and in the pumping frequency, taking into account the Coulomb interaction 
non perturbatively. We find that pumping in a single-dot AB interferometer has a peristaltic but 
phase-coherent character. In a double-dot AB interferometer, we find a pumping mechanism that 
relies purely on quantum-mechanical interference and has no classical counterpart. 

PACS numbers: 73.23.Hk, 85.35.Ds, 72.10.Bg 



I. INTRODUCTION 

Multiply connected mesoscopic structures are ideal to 
study quantum interference in a solid-state system. In 
an Aharonov-Bohm (AB) geometry, the sensitivity of the 
current to the magnetic flux enclosed by the two interfer- 
ing paths can be used as a measure of the quantum coher- 
ence in the system. Embedding quantum dots in the arms 
of the interferometer allows the investigation of the coher- 
ence of transport through a region with strong Coulomb 
interaction. Several experiments 1 - - — have confirmed that 
the visibility of AB oscillations in the current is not com- 
pletely suppressed by the presence of a quantum dot, in- 
dicating that transport though the strongly-interacting 
dot is partially coherent. One mechanism of decoher- 
ence that has been investigated both experimentally 6 and 
theoretically— is spin-flip tunneling. 

The goal of the present paper is to address the issue 
of coherence in adiabatic pumping through systems with 
strong Coulomb interaction. To this end, we consider 
pumping in AB-interferometer devices with a quantum 
dot embedded either in one or in both arms. 

Pumping is a transport mechanism, which exploits 
the periodic time dependence of some parameters of a 
nanoscale conductor to produce a dc current in the ab- 
sence of an applied bias voltage. Its appeal for both 
theorists and experimentalists lies in the possibilities it 
offers to investigate the non-equilibrium induced by the 
explicit time-dependence of a nanoscale system. The 
adiabatic-pumping regime is characterized by the pump- 
ing frequency being smaller than the characteristic time 
scales of the system. Recently, there have been sev- 
eral experiments on pumping in nano systems A lot 
of the theoretical effort has been devoted to systems, 
where the Coulomb interaction can be treated within 
a mean-field approach. 16-25 In this regime, a well es- 
tablished theoretical framework for pumping, based on 
the dynamical scattering approach to mesoscopic trans- 
port, cxistsi 16 ' 26 However, in some nano-scale systems, 
such as few-electron quantum dots, Coulomb interaction 
can become very important, requiring a non-perturbative 



treatment. In the last few years, pumping in strongly 
interacting systems has attracted a lot of theoretical 
interest iSZr— In the present paper, we employ a diagram- 
matic approach to adiabatic pumping in quantum dots, 34 
which relies on a systematic expansion in both the pump- 
ing frequency and the tunnel-coupling strengths. This 
formalism is valid in the weak-tunneling regime but it 
takes into account the on-site Coulomb repulsion in the 
dot non-perturbatively. 

Sometimes the term adiabatic quantum pumping is 
used in the literature to emphasize the role of quantum 
interference in a given pumping mechanism, in contrast 
to a purely classical pump. This motivates the question 
whether such a distinction is well defined. We find that 
this is not always the case. For, example, the pump- 
ing mechanism we discuss in the case of a single-dot AB 
interferometer exhibits both quantum and classical fea- 
tures at the same time. For the case of a double-dot AB 
interferometer, on the other hand, we identify a pumping 
scheme that relies exclusively on quantum-mechanical in- 
terference. 

An important issue when interpreting experimental 
data is to distinguish pumping from rectification4i In 
fact, due to the presence of stray capacitances, undesired 
ac bias voltages may appear across the time-dependent 
conductor, and can give rise to a dc current. 

In this paper, we compare pumping and rectification, 
analyze the different processes contributing to transport 
and discuss to which degree symmetry with respect to 
the magnetic field can be used to distinguish the two 
transport mechanisms. 

The paper is organized as follows. In Sec. [Ill we present 
the model and the technique used to compute the pump- 
ing current. The results for an interferometer with one 
quantum dot embedded in only one of the arms and with 
a quantum dot in both arms are presented in Sees. 1111 Al 
and IIIIB1 respectively. Finally, conclusions are drawn in 

Sec. En 
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FIG. 1: Setup of (a) single-dot and (b) double-dot Aharonov- 
Bohm interferometer. 

II. MODEL AND THEORETICAL APPROACH 
A. Model 

We start by defining the different building blocks of the 
quantum-dot Aharonov-Bohm interferometers depicted 
in Fig. [T] The quantum dots, numbered by the index 
j, are assumed to be in the single-level regime, i.e., they 
can be viewed as Anderson impurities, 

H dottj = Ej n 3<y + U n n n oi . (1) 

Here, rij a — d^ a dj a , with dj a being the creation operator 
for an electron with spin a in quantum dot j. The dot- 
level position is denoted by tj and the onsite Coulomb- 
repulsion energy by U. The inter-dot charging energy 
for a double-dot interferometer is assumed to be negli- 
gible. The two leads are modeled as reservoirs of non- 
interacting electrons, 

Pleads ^ ^ £rkC r k a Crk(7 ) (2) 
rka 

where c\ kty is the creation operator for an electron in lead 
r = L, R in a state labeled by the quantum number k and 
with spin a. 

The tunnel coupling between dot j and the leads is 
modeled by the tunneling Hamiltonian, 

fftunn.j = tr J ( C lkA° + H - C - ) • ( 3 ) 

kor 

We assume the tunnel matrix elements t r j and the den- 
sity of states N r in the lead r to be energy independent 
in the energy window relevant for transport. Tunnel- 
coupling strengths are then defined as T r j = 2ir\t r j\ 2 N r . 
Furthermore, we define Tj = ^2 r T r j. 

An interferometer arm without a quantum dot (refer- 
ence arm) is modeled by a direct tunnel coupling between 



the leads, 

Hrci= fcLo^ + H - c -) » ( 4 ) 

k£R,q£L,<j 

with transmission amplitude t rcl — 27TV^Vl1Vr t. 
In this paper we focus on two different setups: 

(i) A single-dot Aharonov-Bohm interferometer, as 
shown in FigO] (a). In this case, in order to avoid clut- 
tering the notation, we consistently drop the dot index 
(j = 1). The Hamiltonian is simply H — H^ot + H tU nn + 
-f^ieads + Href- The magnetic flux $ threading the in- 
terferometer is included in the phases of the tunneling 
amplitudes. We choose the gauge in which tz,,tji € 5R + 
and argi = <\> = 27r$/<I>o, where $ is the magnetic flux 
and $o is the flux quantum. 

(ii) A double-dot Aharonov-Bohm interferometer, as 
shown in FigJT] (b). In this case, the Hamiltonian reads 
H = Ej=i,2 [ H dot,j + fftunnj] + Pleads- We choose 
the gauge in which — arg^Li = argiru = — argiR2 = 
argt R i = </>/4. 



B. Real-time diagrammatic approach to pumping 

We generalize the real-time diagrammatic approach to 
pumping introduced in Ref. H3 for a single-dot or double- 
dot Aharonov-Bohm interferometer. New ingredients are 
the possibility of direct tunneling from source to drain for 
the reference arm in the single-dot case and the inclusion 
of off-diagonal elements of the reduced density matrix 
which account for coherent superpositions of the elec- 
trons in different quantum dots in the double-dot case. 

The systems under consideration can be decomposed 
in a subsystem with few degrees of freedom comprising 
the quantum dots, and the leads which possess a large 
number of non-interacting degrees of freedom. Since we 
are not interested in the dynamics of the leads' degrees 
of freedom, we can trace them out thus obtaining an ef- 
fective description of the quantum dots' subsystem. The 
Hilbert space of the reduced system is spanned by the 
eigenstates \x) of the Hamiltonian of the isolated dot(s), 
Ylj -Helot, j- The corresponding eigenenergies are denoted 
by E x . For the single-dot natural choice for the 

basis {|x» is {|0>, | t),l \ U) = 44l°»' where the 
different states correspond, respectively, to empty, singly 
occupied by spin up, singly occupied by spin down, and 
doubly occupied dot. The dynamics of the dots' degrees 
of freedom are fully described by the reduced density ma- 
trix p, whose matrix elements arep** = (xi|p|X2)- A di- 
agonal element of the reduced density matrix denotes 
the probability of the dot being in a state x- We in- 
troduce the vector 7r = {p*\, ...,p*™, , •■•) T , (with 
i j), whose first m components are all diagonal ele- 
ments of the reduced density matrix followed by all off- 
diagonal elements. The dynamics of reduced system is 
governed by the generalized master equation (in matrix 
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notation) 

-^Tr(t) = -iE(t)w(i) + /" dt'W(t,t')Tv(t'), (5) 



where the Kernel element W*,£„ (t, t') describes the tran- 
sition from an initial state described by p*.,,, to a final 
state described by p*,. The matrix elements of E(i) are 

given by = <W"W" (fx(*) " 

We are interested in pumping, i.e., in transport due to 
the periodic variation of the system parameters, collec- 
tively denoted by X. The vector 7r(t) as well as the kernel 
W (£,£') depend in a functional way on the pumping pa- 
rameters X(t). To solve Eq. ([5]) we perform an adiabatic 
expansion, 34 i.e., an expansion in powers of the pumping 
frequency f2, which is valid when the pumping frequency 
is much smaller than the response time of the system. For 
this purpose, we first perform a Taylor expansion around 

the final time t of tt(t) = tt(*) + (t - t)^{r) 

integral on the right-hand side of Eq. ((5|). Furthermore, 
we need to perform the adiabatic expansion of the ker- 
nel W(t, t'). In order to do so, we expand the parameters 

around the time t, i.e., X{t) =X(t) + {r- t)^pX{r) 

We write the kernel expansion as 



W(i, *') = W t (i) (t - t') + Wj aJ (t - t') 



(a). 



(6) 



The subscript t denotes the time t around which the adi- 
abatic expansion is performed. The instantaneous part 
[with superscript (i)] is obtained by freezing all parame- 
ters at time t. The adiabatic correction term [with super- 
script (a)] contains only terms which are linear in time 

derivatives of the pumping parameters -4-X(t) . Fi- 

T — t 

nally, we perform an adiabatic expansion of the reduced 
density matrix, 



TT(t) 



» 



(7) 



The instantaneous part can be obtained by solving the 
generalized master equation in the stationary limit 







(-iE(i) 



W 



(8) 



(i) 



1 



together with the normalization condition rnr 
with n = (1, ...1, 0, 0), i.e., the first m components 
of n are 1 and the other components arc 0. In Eq. ©, 
we have introduced the generalized rates as the Laplace 
transform of the Kernel computed at zero frequency: 



W 



(i/a) 



= lim 



/ dt'e-^-^W^it-t'). 



The first adiabatic correction of the generalized master 
Eq. © reads 



dt 



7r« = (-iE(t)+Wr)7r 



w 



( a U(i) 



aw 



(i) 



r (i) 
(9) 



with <9Wf } 



lim j- J dt'e-^*-*')W ( (i) (t-t')- Equs 



tion Q together with the normalization condition 

fa) 

nir t = allows to determine the adiabatic correction 

of the reduced density matrix 7r j '' . 

In the following, we concentrate on the limit of weak 
tunnel couplings. Therefore, we perform a perturbation 
expansion in the tunnel-coupling strength T between dot 
and leads. The fcth order contribution to the reduced 
density matrix is denoted by ■n^ a ' k \ Matching the or- 
ders in Eq. ©, it is easy to see that the expansion of the 

instantaneous term of the reduced density matrix Tv[ hk ^ 
starts in zeroth order (fc = 0), while the adiabatic correc- 
tion ir^'^ starts in minus first order in T {k = — 1). This 
does not invalidate the expansion, since due to the low- 
frequency condition fi -C Y the correction 7r| a ' ^ oc fl/T 
still remains small. In the single-dot AB interferome- 
ter, we also need to consider direct tunneling processes 
between the two leads. This is done by performing a 
in the perturbation expansion in \t \. The order of |t rcf | is in- 
dicated by another superscript I, i.e., the corrections to 
the reduced density matrix are now denoted by ■ K ^/ a,,k,l \ 
The expectation value of the current flowing into lead 
r consists of an instantaneous part and its adiabatic cor- 
rection. The instantaneous part reads 



J»(f) = enW t r ' (i) 7r (i) 



(10) 



with the current rates W r . The latter are calculated sim- 
ilar to W but are weighted with the number of electrons 
transferred to lead r. Without applied bias voltage the 
instantaneous part vanishes. Hence, the pumped current 
is given by the adiabatic correction 



jW(t) = en ( W;'< (a) 7r (i) 



.W r t ' {[) 7Z { t a) 



awr (i) -i 

* dt 



where the superscript r points out that the rates, both 
the instantaneous ones and their adiabatic corrections, 
are current rates and the number of transferred electrons 
needs to be accounted for. 



III. RESULTS 

Using the real-time perturbation theory outlined in 
the previous section, we compute the pumped current 
through an AB interferometer with a quantum dot em- 
bedded either in one or in both arms. We concentrate 
on the limit of weak tunnel coupling and we expand the 
pumped current up to the lowest order that is AB flux de- 
pendent, which is associated with AB interference. Fur- 
thermore, we calculate the dc current driven through the 
AB interferometer by an applied ac bias voltage and rec- 
tified by the time dependence of the instantaneous con- 
ductance of the system. We identify the characteristic 
features for both transport mechanisms. This helps us 



4 



to deepen our understanding of pumping, but more im- 
portantly, to distinguish pumping from rectification in an 
experiment. 



A. Single-dot AB interferometer 

1. Weak adiabatic pumping 

Without applied bias voltage, the instantaneous part 
of the current vanishes. The lowest order of the adiabatic 
correction, 



J, 



(a, 0,0) 
L 



(*) = ~e 



r dt 



>(i,0,0) 



(12) 



is proportional to the time derivative of the average dot 



occupation (n 



.(i,o,o) _ 



2/(6) 



where /(e) is the 



l+/(e)-/(e+t/) ' 

Fermi distribution. Equation (|12[) has very simple inter- 
pretation: if the average occupation changes as a result 
of varying the dot level e, then charge flows into or out 
of the dot. The fraction of the current flowing through 
the left lead is simply given by the ratio 



L.(i 



,1,0) 
X 



W, 



(1,1,0) 



r 



(13) 



where W^l^f^ is the golden-rule rate for a transition 
from x to x' with an electron tunneling through the 
left barrier, while W^/^\^ is the total rate for tunnel- 
ing through the left and right barriers. Notice that in 
the lowest-order perturbation theory considered here, the 
golden-rule rates coincide with the generalized rates in- 
troduced in the previous section. The ratio Eq. (|13l) is 
independent from the initial and final states, x an d x' ■ 
The first flux-dependent correction to the current is 



(a, 0,1) 



(*) 



.VTiTW 

! r 



rcf 



' dt 



v 0,0,0) 



(14) 



It can be interpreted in a similar way as the lowest-order 
contribution with the only difference that now the ratio 
between the tunneling rates is given by 



W L,0,1,1) 

w (i,i,o) 
x'^x 



r 



ref 



Those rates describe the flux-dependent parts of the pro- 
cesses that fill or empty the dot. The flux dependence 
arises due to the interference of the two possible paths 
available: either direct tunneling between dot and lead r 
or the indirect paths that transfers an electron between 
dot and lead r via the other lead f. These rates, that 
change the dot occupation, have quite different proper- 
ties than the rates of same order that describe transfer 
of electrons from one lead to the other.— In particular, 
they are odd functions of the AB flux. 



In the adiabatic regime, to obtain a non-vanishing 
pumped charge at least two parameters (X\ and X2) 
need to be time dependent. We write the parameters 
as Xi(t) — Xi + SXi(t), where Xi is the mean value 
and SXi(t) is the oscillating component. We indicate the 
pumped charge due to the variation of X\ and Xi as 

27r/f2 

Qx 1 ,x 2 'i it can be compute as Qx u x 2 = J dtJ^it). In 



the following, we consider weak pumping and compute 
the pumped charge in bilinear order in 5Xi(t). Since the 
current is proportional to the time derivative of the aver- 
age dot occupation, see Eqs. ([i"2"j) and ([T4"]) . one pumping 
parameter needs to be the dot level e. Choosing IY to be 
the second pumping parameter we obtain for the pumped 
charge in zeroth order in V and in zeroth and first order 
in |i rof | 



Qr 



(a, 0,0) 



r R d 
' e f^ vr ^di 



■ 0,0,0) 



(15) 



n (a,0,l) 



where the prefactor 




>— 

de 



,(i,0,0) 



(16) 



2tt/Q 



—ST r dt 



is the area of the pumping cycle in parameter space. The 
charge for pumping with (Tr,, e) is simply obtained from 
Eq. (fT5"|) swapping IY and Tr . 

The dependence of the pumped charge on the aver- 
age level position in zeroth order in T and up to first 
order in \t Toi \, see Eqs. (TT51) and (IT^l) . is simply given 
by 3! (n) (i,0 ' 0) . The latter is plotted in Fig. The 
pumped charge is even around s = —{7/2. The direction 
of the pumped current is independent of the average dot 
level. The two peaks in Fig. [5] are associated to the 
transitions between singly occupied and empty dot and 
between doubly occupied and singly occupied dot. 

Is the mechanism of pumping in our example of classi- 
cal or of quantum nature? The interpretation of Eqs. (|12[) 
and (I14j) is consistent with a picture of a peristaltic pump: 
a variation of the gate voltage pushes the electrons off the 
dot, which generates a current flow to the leads, where 
the larger part flows through the contact that is more 
open, as characterized by the relative coupling strengths 
to the left and right leads. If this ratio is changed when 
the gate voltage sucks electrons into the dot later, then 
a net charge has been transferred from one lead to the 
other. This mechanism does not rely on the formation of 
a superposition of quantum states. We refer to it as peri- 
staltic, which is usually associated with classical pump- 
ing. On the other hand, the underlying tunnel processes 
that transfer the dot electrons from or to the leads are 
phase coherent, as signaled by the dependence on the AB 
phase. They are a consequence of quantum-mechanical 
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FIG. 2: The dependence of pumped charge in zeroth order 
in r and up to first oder in |i ref | on the average dot level 



e is given by the derivative of the average dot occupation 
21 (n)^'°'°\ which is plotted here for various Coulomb inter- 
action strengths U . The temperature is k^T — 2T. 



interference of the two possible paths between dot and a 
given lead. Therefore, we conclude that a clear distinc- 
tion between quantum and classical pumping is mean- 
ingless for our example. To emphasize the coexistence 
of both classical and quantum features, we describe the 
mechanism considered here as phase-coherent peristaltic 
pumping. 



2. Comparison with rectification 

Adiabatic pumping may be obscured by rectification. 
A time-dependent gate voltage may not only change the 
level position in the quantum dot but, due to a parasitic 
capacitive coupling to the leads, give rise to an effec- 
tive (in-phase) ac bias voltage. This ac bias voltage can, 
in turn, yield a dc current component due to the time 
dependence of the dot-level position. In Ref. 0, sym- 
metry with respect to magnetic held has been used to 
discriminate pumping from rectification. Assuming that 
the lever arms between gates and reservoirs are small, 
one can neglect rectification contributions quadratic in 
V(t) but to zeroth order in time variation of the sys- 
tem parameters of the pumping region (Se in our case). 
This is because the effect of the gate- voltage modulation 
on the dot-level position dominates over the ac voltage 
due to the parasitic stray capacitance. In this limit, the 
charge transferred in one period by rectification can be 
computed as Q lcc ,x = 

J^ /n G h (t)V(t)dt, where G L (t) 
is the instantaneous linear conductance and V(t) is the 
undesired oscillating bias voltage. Due to Onsager re- 
lations the linear conductance, and, therefore, also the 
rectification contribution to the transferred charge, is an 
even function of the magnetic held. This reasoning, how- 
ever, is no longer valid when contributions to the recti- 
fied current that are non-linear in the parasitic ac bias 
voltage have to be taken into account. In fact, mag- 



netic held symmetries for different transport regimes have 
been extensively investigated experimentally 4 ^— as well 
as theoretically ! 48 ' 49 In nonlinear response it has been 
measured that Coulomb interaction may yield an odd 
part also in rectification.— The ratio between odd and 
even parts strongly depends on the bias mode and the 
frequency. Especially, in the adiabatic regime the odd 
part is in general not negligible— 

In the following, we choose e as time-dependent pa- 
rameter and we compute the rectified charge in linear 
order in Ss(t) and V(t). In lowest non-vanishing order in 
F, the charge transferred by rectification reads 



Q 



(M,0) = 2 
rcc.e ° 



r L r r 



-»?r. 



d_ 

els 



r 

(1 - f(e + U)) fj{e) + f(g)*f(e + U) 



l + f(e)-f(e + U) 



(17) 



with n re 



J^ /n dtSe(t)V{t). The rectified charge 



Eq. (T7]) is odd around e = -U/2 (see Fig©. 
The first flux-dependent correction reads 



rcc.e 



with 

/^(i.broad) _ 



2 2 vTLTRuref 



cos 0— (n 
ale 



(i, broad) 



(18) 



= 2- (n 



,(i,0,0) 



2tt de 



Re 



,0,0,0) 



r d 

2nde' 



-Re 



1 

2 " 
ip 
2n 



IB 
(e + U) 



where W is the digamma function. Unlike the lowest non- 
vanishing order the first flux-dependent correction of the 
charge is even around e — (see Fig [4]). 
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FIG. 3: Charge transferred by rectification in lowest non- 
vanishing order , Qrec.'e , in units of Qq — e 2 F pj L n rcc>£ as 
a function of the average dot level e for various Coulomb- 
interaction strengths U. The time- varying parameter is e and 
the temperature is feT = 2f. 

In the linear-response limit, the transferred charge due 
to rectiheation is even with respect to magnetic flux cj). 
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FIG. 4: First flux-dependent correction to the charge 
transferred by rectification, Qrec.'e ! in units of Qo = 
„2 \Z r Lr R I 



ef I 



i _• ■ |?7rec,e cos cf> as a function of the average dot level 
e for various Coulomb-interactions strengths U . The time- 
varying parameter is e and the temperature is k^T = 2F. 



On the contrary, the flux-dependent part of the pumped 
current is odd. The different symmetry with magnetic 
field is related to the different processes which contribute 
to pumping and to rectification. For rectification, the 
processes that contribute to transport describe charge 
transfer from one lead to the other. The flux- dependent 
parts are associated with the interference of cotunneling 
through the dot and direct tunneling through the ref- 
erence arm. This is different for pumping, as discussed 
in the previous section. There, the relevant processes 
change the dot occupation and are resonant with inter- 
mediate state on one of the leads. They are associated to 
interference between direct tunneling from the dot to a 
given lead and cotunneling from the dot to this lead via 
the other lead. 

Furthermore, the symmetry with respect to e is differ- 
ent for pumping and rectification. The pumped charge 
is even about e = —U/2, while for rectification there is a 
different symmetry in different orders of the perturbation 
expansion. 



B. Double-dot AB interferometer 

For the double-dot interferometer, we consider two dif- 
ferent limits regarding the Coulomb-interaction strength: 
(i) fully non-interacting case and (ii) infinite intra-dot 
interaction, which forbids double occupation of a single 
dot, and negligible inter-dot interaction. 



1. Weak adiabatic pumping 

We consider the tunneling barriers between dot and 
lead to be the same for both dots: IYi = Tl2 = IY and 
Tri = TR2 = Tr. We assume the difference Ae = e\ — e-i 
between the dot level of the upper and of the lower dots 



to be of the same order as T. The average level is de- 
fined as e = (ei + £2)/2. We calculate the current in 
first order in T. This order of perturbation theory is 
already flux dependent. In the non interacting case, we 
can consider spinless electrons and take into account spin 
degeneracy by multiplying the current by a factor of 2. 
The dot Hilbert space for spinless electrons is spanned by 
the states {|0), |1), |2), |12)}, corresponding, respectively, 
to both dots being empty, only upper dot occupied, only 
lower dot occupied, and both dots occupied. On the other 
hand, for infinite intra-dot interaction the dot Hilbert 
space has dimension 9 and it is spanned by the states 
{\0),\ja),\la2a')}, with j = 1,2 and a, a 1 =t,|. These 
states correspond, respectively, to both dots being empty, 
dot j occupied with spin a and both dots occupied with 
spin a in dot 1 and spin a' in dot 2. For the interact- 
ing system, we introduce the abbreviation = pj* . In 
both cases addressed here, the pumped current is writ- 
ten conveniently as a function of the isospin's expectation 
value 



1=1 h 



1 i 2 
1 / P2 + Pi 

= 5 i(P?-P&) 

v\-v\ 



Notice that if the isospin lies in the xy plane, it indicates 
that the system is in a superposition of two states: one 
with only an electron in the first dot and a second with 
only an electron in the second dot (both states with the 
same spin). 

Computing the adiabatic correction to the reduced 
density matrix and inserting it in Eq. in isospin no- 
tation, we obtain for noninteracting electrons 

(a0) Fl 2r R .sinf (rsinf + Aecosf)+Ae 2 

J L(^o)- 4e r 4r L r R sin 2 |+A £ 2 



r 

df{e) de 
de dt 



(19) 



In the case of infinite interaction within one dot and van- 
ishing interaction between the dots the expression for 
the current is quite long. For symmetric tunnel-coupling 
strengths (Fl = Tr_) it simplifies to 



J, 



(a,0) 

L([7-+oo) 



-2e 



df(s) de 
de dt 



T 2 sin 2 | [1 + /(e)] 3 + Ae 2 [l + /(e)] + TAe cos § sin f 



[l + /(e)]3{(r 2 sin|[l + /(e)] 2 + Ae 2 } 



(20) 



Now we consider weak pumping with the parameters 
e(t) = e + 5e(t) and Ae(t) = Ae + SAe(t). The area of 
the cycle is 



2tt/Q 



dSe 
dt 



SAedt. 
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For noninteracting electrons the pumped charge per pe- 
riod reads 



Q 



(a,o) o r LrR . (j)df(e) 
- 8e?7A e , e „ sin 



Ae,e(!7=0) '" e > £ p 2 rfe 

cos f ( Ae" 2 - 4r L r R sin 2 |) - 2(r L - r R ) sin I Ai 



4r L r R sin 2 ?+ Ae' 



(21) 



For an infinite Coulomb interaction we give, again, the 
expression for symmetric tunnel coupling, 



Q (.,0) 



Ae,e(f7->oo) 



2er? 



Ae,e-l 



ds 

2 



cos I sin |{ Ae 2 - T 2 sin 2 § [1 + /(e)] 2 } 



[l + /(e)] 3 {r2sin 2 |[l + /(e)] 2 + Ae 2 } 

(22) 

The pumped charge for Tl ^ r R as a function of e and 
Ae is shown in Figs.jSJa) and 5(b). We find a sign change 
in the pumped charge which in the noninteracting case 
only depends on Ae but for an infinite interaction also 
depends on e . Equation [2D suggests an even symme- 
try concerning Ae. Figure O:) shows that this symme- 
try is not general but only valid for symmetric tunneling 
barriers. As a function of e, the pumped charge is even 
only in the noninteracting but not in the interacting case, 
Fig. Ed). 

The fact that we find a nonvanishing pumped charge at 
all is not self-evident. The two pumping parameters are 
associated with the different arms of the interferometer. 
This suggests that pumping relies on coherent superpo- 
sition of states localized in the different arms described 
by the isospin components I x and I y . Therefore, one can 
view pumping in this case as fully quantum mechanical. 



2. Comparison with rectification 

Similarly to the case of a single-dot interferometer, 
we consider rectification in the linear-response regime, 
in which the linear conductance, and, therefore, also the 
transferred charge, is an even function of the magnetic 
flux. The linear conductance for vanishing interaction 
reads 



G 



(i) 

L(E/=0) 



lc2 r L r R ,d/ Ae 2 + r L r R sin^ 

r de Ae 2 + r L r R sin 2 | ' 1 ' 



while for infinite intra-dot interaction it is 

2 r L r R i df 



G 



L(C/^oo) - r i + f( £ -)de 
r L r R (l - cos 0) [1 + cos 4> + 2/(e) (2 + /(e))] + Ae 2 



Q (a ' 0) /Q 



2r L r R (l - cos0) (1 + /(e)) 2 + Ae 2 



(24) 




Q (a ' 0) /Q 




O 



0.08 






0.07 




0.06 




0.05 


Q 




Q> 


0.04 


O oj 




— < 


0.03 








0.02 




0.01 








FIG. 5: Pumped charge Q^fl in units of Qo — erj(e, <5e)/T 2 
for T L = 0.8r, T R = 0.2r, <j> '= tt/2, and k B T = 2F. (a) and 
(b) show a density plot where Q&fl is a function of e and Ae 
for (a) vanishing and (b) infinite Coulomb interaction. Cuts 
through (b) are shown in (c) and (d). In (c) Q^f\ is plotted 
as a function of Ae for different e. In (d) Q^f\ is plotted as 
a function of e for different Ae. 
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® 0.01 



-0.01 



- \ / 




\ / - 






pumping 

rectification | 



FIG. 6: Pumped charge Q^fl in units of Qo — erj(e, fe)/P 2 

and rectified charge Qx<x}e in units of Qo = e 2 ri ICCy£ /T as a 
function of for U = oo, T L = 0.8r, T R = 0.2r, e = 0, 
Ae = 0.5r, and k B T = 2T. 



ther in one or in both arms, by means of a diagram- 
matic real-time approach to pumping. In the single-dot 
AB interferometer, we have found that adiabatic pump- 
ing has a peristaltic character. Nonetheless, it is clearly 
phase coherent as indicated by the flux dependence of the 
pumped current. On the other hand, in a double-dot AB- 
interferometer adiabatic pumping with the levels of the 
two dots is a pure quantum-mechanical transport mech- 
anism, since it relies on the system being in a coherent 
superposition of eigenstates of the dots in the upper and 
lower arms. This pumping mechanism has no classical 
counterpart. Finally, we found that the symmetry of the 
pumped charge with respect to the magnetic flux may 
help to distinguish pumping from rectification, at least 
in the linear-response regime. 



The pumping, on the other hand, has no definite sym- 
metry with respect to magnetic field (Fig. [5]) unless a 
symmetric choice of the tunnel-coupling strengths is as- 
sumed. Furthermore, we remark that the pumped charge 
vanishes for zero flux. 



IV. CONCLUSIONS 

We have investigated adiabatic pumping through an 
AB interferometer with a quantum dot embedded ei- 
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